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Be s chr e ibung 

Rechnergestiitztes Modellierverf ahren fiir das Verhalten eines 
Stahlvo lumens mit einer Volumenoberf lache 

Die vorliegende Erfindung betrif ft ein rechnergestiitztes Mo- 
del lierverfahr en fiir das Verhalten eines Stahlvolumens mit 
einer Volumenoberf lache, 

- wobei ein Rechner anhand eines moment anen Anf angszustands 
des Stahlvolumens und mindestens einer iiber die Volumen- 
oberflache auf das Stahlvolumen einwirkenden momentanen 
Einf lussgroSe durch Losen einer Warmeleitungsgleichung und 
einer Phasen\jmwandlungsgleichung einen Folgezustand des 
Stahlvolumens ermi ttelt , 

- wobei die mindestens eine Einf lussgroSe fiir eine Anzahl von 
Flachenelementen der Volxamenoberf lache jeweils einen loka- 
len Einfluss umfasst und die lokalen Einflusse iiber das je- 
weilige Flachenelement auf das Stahlvolumen einwirken, 

- wobei der Anf angszustand und der Folgezustand fiir eine An- 
zahl von Volumenelementen des Stahlvolumens jeweils lokale 
Anteile von modellierten Phasen des Stahls und eine einen 
lokalen Energieinhalt des Stahls beschreibende GroSe iimfas- 
sen, 

- wobei die modellierten Phasen des Stahls Austenit und eine 
erste weitere Phase umfassen, in die Austenit ijinwandelbar 
ist und die in Austenit umwandelbar ist. 

Ein derartiges Modellierverf ahren ist z. B. aus der DE~A-101 
29 565 bekannt. Bei dieser Schrift wird insbesondere erstmals 
versucht, die Fouriersche Warmeleitungsgleichung selbst und 
nicht eine unkorrekte Abwandlung dieser Warmeleitungsglei- 
chung zu losen, um das thermodynamische Verhalten eines 
Stahlbandes korrekt zu beschreiben. Diese Schrift wird hier- 
mit durch Verweis auf sie in den Of f enbarungsgehalt der vor- 
liegenden Erfindung mit aufgenommen. 
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In der alteren, zum Anmeldezeitpunkt der vorliegenden Erfin- 
dung nicht vorverof f entlichten deutschen Patentanmeldung 102 
51 716.9 ist ebenfalls ein derartiges Modellierverf ahren be- 
schrieben. Bei diesem Modellierverf ahren wird versucht, die 
5 Phasenumwandlung des Stahls anhand der Gibbsschen freien 

Enthalpien des Stahls zu model lier en. Auch diese Schrift wird 
hiermit durch Verweis auf sie in den Of f enbarungsgehalt der 
vorliegenden Erfindung mit auf genommen. 

10 Weiterhin ist ein derartiges Modellierverf ahren z. B. aus dem 
Aufsatz ^^Numerische Simulation der Warmeleitung in Stahlble- 
chen - Mathematik hilft bei der Steuerung von Kiihlstrecken'' 
von W. Borchers et al., verof f entlicht im Unikurier der 
Friedrich-Alexander-Universitat Erlangen-Niirnberg, Heft 102 , 
15 Oktober 2001, 27. Jahrgang, bekannt. 

SchlieSlich sind noch die traditionellen Ansatze gemafi der 
Scheilschen Kegel, nach Johns on-Mehl-Avr ami und Brimacombe 
bekannt . 

20 

Bereits die exakte Modellierung des zeitlichen Temperaturver- 
laufs von Stahl bei der Abkiihlung, insbesondere von Stahlban- 
dern, ist entscheidend fiir die Steuerxang der erf orderlichen 
Wasser- bzw. Kiihlmittelmengen einer Kiihlstrecke fur Stahl. 
Denn die beim Abkiihlen auf tretende Stahlumwandlung beein- 
flusst entscheidend das thermische Verhalten des Stahls bei 
der Abkiihlung. Auch warden durch den Abkuhlprozess wesentli- 
che Materialeigenschaf ten des Stahls beeinflusst. Da die Ab- 
kiihlung aber nicht im thermodynamischen Gleichgewicht er- 
3 0 folgt, ist es nicht moglich, die Umwandlung einfach durch ge- 
eignete Anpassung der Warmekapazitat zu beschreiben. Daher 
ist zusatzlich eine exakte Modellierung der Phasenumwandlung 
des Stahls erf orderlich, um die Kiihlstrecke korrekt steuern 
zu konnen. 



35 



Die traditionellen Ansatze des Standes der Technik arbeiten 
in der Praxis nicht in alien Fallen fehlerfrei. Insbesondere 
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weisen sie eine Reihe von systematischen Nachteilen auf . Zum 
Ersten muss jedes Material separat parametriert werden. In- 
terpolationen zwischen verschiedenen Materialien sind nicht 
Oder zumindest nur eingeschrankt moglich. Zum Zweiten werden 
5 bei den traditionellen Verfahren des Standes der Technik nur 
zwei Phasen betrachtet. Eine Erweiterung auf mehr als zwei 
Phasen ist systembedingt nicht moglich. Zum Dritten liefern 
die Verfahren des traditionellen Standes der Technik nur bei 
einer vollstandigen Umwandlung des betrachteten Me tails eine 
10 gute Ubereinstimmung zwischen Modell und Realitat. Zum Vier- 
^ ten liefern die Verfahren des traditionellen Standes der 
M Technik keine Information iiber die bei der Phasenumwandlxang 
^ freigesetzte Umwandlungswarme . Die Kenntnis der Umwandlungs- 
warme ist fur ein korrektes Losen der Warmeleitungsgleichung 
15 aber unabdingbar. 



Die Vorgehensweisen gemaS der DE-A-101 29 565 und dem Fach- 
auf satz ^'Niomerische Simulation . . . " stellen dem gegenuber be- 
reits einen erheblichen Fortschritt dar, da sie zumindest die 

2 0 Warmeleitung vollig korrekt beschreiben. Die altere deutsche 

Patentanmeldung verbessert dariiber hinaus die Modellierung 
der Phasenumwandlung . Insbesondere liefert sie bereits die 
bei der Phasenumwandlung auftretende Umwandlungswarme. Auch 
^ diese Verfahren sind aber noch verbesserungsf ahig. 

Die Aufgabe der vorliegenden Erfindung besteht darin, ein Mo- 
dellierverf ahren fur ein Metall zu schaffen, das bessere Mo- 
de 11 ergebnisse liefert. 

3 0 Die Aufgabe wird dadurch gelost, 

- dass der Anf angszustand und der Folgezustand fiir mindestens 
eines der Volumenelemente auch eine lokale Konzentrations- 
verteilung mindestens eines im Stahl beweglichen Legie- 
r ungs e 1 emen t s umf as s en , 
3 5 - dass im Rahmen der Umwandlungsgleichung fur das mindestens 
eine Volumenelement ermittelt wird, welche Konzentrationen 
des mindestens einen beweglichen Legierungselements beid- 
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seits einer ersten Phasengrenze zwischen Austenit und der 
ersten weiteren Phase vorliegen./ 

- dass durch. Losen eines ersten Stefan-Problems ermittelt 
wird, ob und wie sich die Konzentrationsverteilung des min- 
xes tens einen beweglichen Legierungselements im austeniti- 
schen Bereich des betrachteten Volumenelements andert und 
ob und Tim welches AusmaS sich die erste Phasengrenze da- 
durch verschiebt, und 

- dass die lokalen Anteile der Phasen anhand einer durch das 
AusmaS der Verschiebung der ersten Phasengrenze bestimmten 
Lage der ersten Phasengrenze ermittelt werden. 

Das in Stahl bewegliche Legierungselement ist in der Kegel 
Kohlenstoff . Alternativ oder zusatzlich ist es aber auch mog- 
lich, dass das Legierungselement Stickstoff ist. Es kann auch 
ein anderes Legierungselement sein, das in Stahl bevorzugt an 
Zwischengitterplatzen angeordnet ist. 

Bereits dieses Modellierverf ahren ermoglicht gegeniiber dem 
bisher bekannten Stand der Technik einen deutlichen Fort- 
schritt. Denn mit dieser Vorgehensweise kann - je nach dem, 
ob die weitere Phase Ferrit oder Zementit ist - das Umwand- 
lungsverhalten von Austenit in Ferrit bzw. von Austenit in 
Zementit und umgekehrt sehr realitatsgetreu modelliert wer- 
den. 

Vorzugsweise \amfassen die modellierten Phasen des Stahls auch 
eine zweite weitere Phase, in die Austenit umwandelbar ist 
und die in Austenit umwandelbar ist. Denn dann ist es insbe- 
sondere moglich, 

- dass fur das betrachtete Volumenelement im Rahmen der Um- 
wandlungsgleichung auch ermittelt wird, welche Konzentrati- 
onen des mindestens einen beweglichen Legierungselements 
beidseits einer zweiten Phasengrenze zwischen Austenit und 
der zweiten weiteren Phase vorliegen, 

- dass durch zusatzliches Losen eines zweiten Stefan- Problems 
ermittelt wird, ob und wie sich die Konzentrationsvertei- 
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lung des mindestens einen beweglichen Legierungselements im 
austenitischen Bereich des betrachteten Volumenelements an- 
dert und ob und urn welches AusmaS sich die zweite Phasen- 
grenze dadurch verschiebt, 

- dass die Stef an-Probleme miteinander gekoppelt sind, 

- dass den Phasengrenzen FlachenmaEe zugeordnet sind, 

- dass ein Anteil des der zweiten Phasengrenze zugeordneten 
FlachenmaSes an der Suxnme der Flacheranafie bestiinmt wird und 

- dass die lokalen Anteile der Phasen aucii vom Anteil des der 
zweiten Phasengrenze zugeordneten FlachenmaSes an der Summe 
der FlachenmaSe abhangen, 

Mit dieser Vorgehensweise ist nicht nur die Umwandlung zwi- 
schen Austenit auf der einen Seite und Ferrit oder Zementit 
auf der anderen Seite sehr realitatsgetreu beschreibbar , son- 
dern insbesondere auch die Umwandlung von Austenit in Perlit 
und umgekehrt. Es ist somit das gesamte wesentliche Umwand- 
lungsverhalten von Stahl in festem Zustand korrekt beschreib- 
bar. Ferner ist mit dieser Vorgehensweise ermittelbar, ob 
Perlit gebildet wird oder nicht. 

Der Anteil des der zweiten Phasengrenze zugeordneten Flachen- 
maSes an der Summe der FlachenmalSe kann der art bestimmt wer- 
den, dass die Phasengrenzen stets nebeneinander angeordnet 
bleiben. Alternativ ist es aber auch moglich, diesen Anteil 
derart nachzufuhren, dass die Phasengrenzen aufeinander zu 
streben. Anhand dieses Anteils kann dann auch abgeleitet wer- 
den, ob Austenit nur in die erste weitere Phase, nur in die 
zweite weitere Phase oder sowohl in die erste als auch in die 
zweite weitere Phase umgewandelt wird. 

Prinzipiell ist es moglich, die Ausbildung der ersten Phasen 
grenze dreidimensional zu ermitteln. Dies f order t aber einen 
erheblichen Rechenauf wand . Vorzugsweise wird daher die Rech- 
nung eindimensional vorgenommen. Dies ist einfacher moglich, 
wenn das Volumenelement ein raumfiillender aristotelescher 
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Korper, insbesondere ein Quader, ist. Es ist also vorzugs- 
weise so, 

- dass das betrachtete Volumenelement quaderformig ausgebil- 
det ist und drei Grundabmessungen aufweist, 

- dass die erste Phasengrenze als Rechteck mit einer ersten 
Langsseite und einer ersten Querseite ausgebildet ist und 

- dass die erste Langsseite mit einer ersten der Grundabmes- 
sungen korrespondiert, die erste Querseite parallel zu ei- 
ner zweiten der Grundabmessungen verlauf t und Verschiebun- 
gen der ersten Phasengrenze parallel zu der dritten der 
Grundabmessungen erfolgen* 

Bei Beriicksichtigung auch der zweiten weiteren Phase ist es 
dann so , 

- dass die zweite Phasengrenze als Rechteck mit einer zweiten 
Langsseite und einer zweiten Querseite ausgebildet ist und 

- dass die zweite Langsseite mit der ersten der Grundabmes- 
sungen korrespondiert, die zweite Querseite parallel zur 
zweiten der Grundabmessungen verlauf t und Verschiebungen 
der zweiten Phasengrenze parallel zur dritten der Grundab- 
messungen er f olgen • 

Wenn die Summe der Querseiten der Phasengrenzen ungefahr 
gleich dem 1,5- bis 3-fachen eines kritischen Lamellenabstan- 
des ist, bei dem eine Energiebilanz , die einerseits mit dem 
Verschieben der Phasengrenzen korrespondierende Phasenumwand- 
lungen des Stahls und andererseits mit dem Verschieben der 
Phasengrenzen korrespondierende Anderungen der Flache einer 
Grenzschicht zwischen der ersten und der zweiten weiteren 
Phase berucksichtigt, ausgeglichen ist, kann mittels des Mo- 
dells auch der Lamellenabstand des Perlits ermittelt werden. 
Der Lamellenabstand (das heiSt die Summe der Querseiten der 
Phasengrenzen) kann dabei insbesondere in etwa gleich dem 
Doppelten des kritischen Lamellenabstands sein. 

Trotz der nur eindimensionalen Rechnung liefert das erfin- 
dungsgemaSe Modellierverf ahren sehr realitatsgetreue Ergeb- 
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nisse/ wenn der Anteil an Austenit anhand einer nicht linea- 
ren Funktion der Lage der Phasengrenze bzw, der Phasengrenzen 
ermittelt wird. 

5 Wenn die Entiittlung, in welchen Konzentrationen das mindes- 
tens eine bewegliche Legierungseleinent beidseits der ersten 
Phasengrenze bzw. beidseits der ersten und beidseits der 
zweiten Phasengrenze vorliegt, anhand der Gibbsschen freien 
Enthalpien der Phasen erfolgt, ist die Ermit tiling der Kon- 
10 zentrationen des mindestens einen beweglichen Legierungsele- 
ments an den Phasengrenzen besonders einfach. 




Es ist moglich, stets die gekoppelten Stef an-Probleme zu 16- 
sen und daraus abzuleiten, welche Phasen in welchem Umfang 
15 gebildet werden. In manchen Fallen ist dies sogar unumgang- 

lich. Manchmal ist es aber moglich, vor dem Losen der Phasen- 
umwandlTingsgleichijng anhand der im Anf angszustand bereits 
vorhandenen Phasen und anhand der Gibbsschen freien Enthal- 
pien vorab zu ermitteln, ob sowohl Austenit als auch die 
20 erste weitere Phase vorliegen bzw. ob zusatzlich zu Austenit 
und der ersten weiteren Phase auch die zweite weitere Phase 
vorliegt . 

Im Einzelfall kann die Anzahl von Volumenelementen klein 
sein. Im Extremfall ist sie gleich Eins . In aller Regel um- 
fasst das Stahlvolumen aber eine Vielzahl von Voliimenelemen- 
ten. Wenn das Stefan-Problem bzw. die Stef an-Probleme nur fur 
einen Teil der Voliimenelemente gelost wird bzw. werden und 
die lokalen Anteile der Phasen der anderen Volumenelemente 
3 0 anhand der lokalen Anteile der Phasen des Teils der Volumen- 
elemente ermittelt werden, kann daher die zum Modellieren des 
Verhaltens des Stahlvolumens erf orderliche Rechenleistung 
deutlich reduziert werden, ohne die Aussagekraft der Modell- 
rechnung wesentlich zu verschlechtern. Die Warmeleitungsglei- 
35 Chung hingegen wird in aller Regel individuell ftir jedes Vo- 
lumenelement gelost . 
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Das Modellierverf ahren kann alternativ online und in Echtzeit 
Oder aber offline ausgefiihrt werden. 

Es ist z. B. moglich, 
5 - dass dem Rechner ein Erstzustand und z-umindest eine ge- 
wiinschte Endgrofie vorgegeben werden, 

- dass das obige Modellierverf ahren iterativ angewendet wird, 

- dass der Anf angszustand der ersten Iteration dem Erstzu- 
stand und der Anf angszustand jeder weiteren Iteration dem 

10 unmittelbar zuvor ermittelten Folgezustand entspricht und 




- dass anhand des nach einer letzten Iteration ermittelten 
Folgezustands eine erwartete EndgroSe ermittelt und mit der 
gewiinschten EndgroSe verglichen wird. 



15 In diesem Fall kann das erf indungsgemaSe Modellierverf ahren 
alternativ online und in Echtzeit oder offline ausgefiihrt 
werden . 

Insbesondere bei einer of f line-Ausfiihrung ist es dabei mog- 
20 lich, dass die Einf lussgroSen der Iterationen in ihrer Ge- 

samtheit einem Einf lussgrofienverlauf entsprechen, der Rechner 
auf Grund des Vergleichs der erwarteten EndgroSe mit der ge- 
wiinschten EndgroSe den Einf lussgroSenverlauf variiert und 
^ ausgehend vom Erstzustand das obige Modellierverf ahren erneut 
H. ausfiihrt, bis zumindest die erwartete EndgroSe mit der ge- 
wiinschten EndgroSe korrespondiert . 

Es ist aber auch moglich, dass der Rechner anhand einer aus 
dem Anf angszustand ermittelten AnfangsgroSe und einer ge- 

30 wiinschten FolgegroSe die Einf lussgroSe ermittelt und eine Be- 
einf lussxingseinrichtung derart ansteuert, dass das Stahlvolu- 
men entsprechend der ermittelten Einf lussgroSe beeinflusst 
wird. In diesem Fall muss das Modellierverf ahren online und 
in Echtzeit ausgefiihrt werden. Das Ansteuern der Beeinflus- 

35 sungseinrichtung kann dabei alternativ sofort oder in der 
darauf f olgenden Iteration erfolgen. 
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Weitere Vorteile und Einzelheiten ergeben sich aus der nach- 
folgenden Beschreibung eines Ausfiihrungsbeispiels in Verbin- 
dxmg mit den Zeichnungen* Dabei zeigen in Prinzipdarstellung 
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FIG 1 
FIG 2 
FIG 3 
FIG 4 
FIG 5 
FIG 6 



FIG 7 



FIG 8 
FIG 9 



FIG 10 



FIG 11 
FIG 12 
FIG 13 
FIG 14 
FIG 15 
FIG 16 

FIG 17 

FIG 18 
FIG 19 



eine Kxihlstrecke fiir ein Stahlband, 
das Stahlband von FIG 1 im Detail, 

ein Voliamenelement in perspektivischer Darstellung, 

ein Ablaufdiagramin, 

ein weiteres Ablauf diagraxnm, 

eine Warmeleitungsgleichung und eine Phasenumwand- 
lung s g 1 e i chung , 

eine weitere Warmeleitungsgleichung und eine Phasen- 

umwandlungsgleichung , 

ein weiteres Ablauf diagramm, 

die Gibbsschen freien Enthalpien der Phasen des 
Stahls als Funktion des Kohlenstof f anteils bei einer 
ersten Tempera tur, 

die Gibbsschen freien Enthalpien als Funktion des 
Kohlenstof f anteils des Stahls bei einer zweiten Tem- 
pera tur , 

ein Volumenelement in perspektivischer Darstellung, 
das Volumenelement von FIG 11 in der Seitenansicht , 
ein Volumenelement analog dem von FIG 12, 
ein Volumenelement in perspektivischer Darstellung, 
das Voliamenelement von FIG 14 in der Draufsicht, 
das Volumenelement von FIG 15 im Schnitt gemaS der 
Linie XVI - XVI von FIG 15, 

das Volumenelement von FIG 15 im Schnitt gemaS der 
Linie XVII - XVII von FIG 15, 
ein weiteres Ablauf diagramm und 
ein Detail von FIG 15 . 



35 



Gemafi FIG 1 soil z. B. ein Stahlband 1 gekiihlt werden, damit 
es vorgegebene Endeigenschaf ten annimmt. Hierzu wird das 
Stahlband 1 z. B. in einer Kxihlstrecke 2 mit einem Kiihlmedium 
3 - in der Kegel Wasser - beauf schlagt • 
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Das Stahlband 1 weist eine Bandbreite b und eine Banddicke d 
auf . Es weist ferner eine — prinzipiell beliebige - Bandlange 
auf . Es durchlauft die Kiihlstrecke 2 mit einer Bandgeschwin- 
digkeit v. 

5 

Die Bandgeschwindigkeit v wird kontinuierlich erfasst und ei- 
nem Rechner 4 zugefiihrt, der die Kiihlstrecke 2 steuert. Der 
Rechner 4 ist somit in der Lage, in an sich bekannter Weise 
eine Wegverf olgung des Stahlbandes 1 zu realisieren. Dem 
10 Rectmer 4 ist also bekannt, auf welche Bereiche des Stahlban- 
des 1 eingewirkt wird, wenn das Stahlband 1 mittels Aufbrin- 
gungseinrichtungen 5 der Kiihlstrecke 2 mit dem Kiihlmedium 3 
beaufschlagt wird. 

15 Der Rechner 4 ist mit einem Computerprogramm 5 programmiert , 

das dem Rechner 4 uber einen Datentrager 7 , z . B . eine CD-ROM 
1 , zugefuhrt wird. Auf dem Datentrager 7 ist das Computerpro- 
gramm 6 in (ausschlieBlich) maschinenlesbarer Form gespei- 
chert. Das Computerprogramm 5 wird vom Rechner 4 iibernommen 

2 0 und in einem Massenspeicher 8, z. B. einer Festplatte 8, des 
Rechners 4 hinterlegt. Dadurch ist der Rechner 4 in der Lage, 
bei Aufruf des Computerprogramms 6 ein Modellierverf ahren fiir 
das Stahlband 1 bzw. fiir dessen einzelne Bereiche (= das 
Stahlvolumen 1) auszufiihren, das nachfolgend naher beschrie- 
ben wird. 

GemaS FIG 2 wird hierzu das Stahlvolumen 1 innerhalb des 
Rechners 4 in Volumenelemente 9 zerlegt. Soweit ein Voliimen- 
element 9 dabei nicht allseitig von anderen Volumenelementen 

3 0 9 umgeben ist, sind dem jeweiligen Volumenelement 9 ein oder 

zwei Flachenelemente 10 zugeordnet. Die Flachenelemente 10 
bilden in ihrer Gesamtheit eine Volumenoberf lache des Stahl- 
bandes 1 bzw. des betrachteten Bereichs des Stahlbandes 1. 

3 5 Der Rechner 4 realisiert auf Grund der Programmierung mit dem 
Computerprogramm 6 insbesondere ein Modell, in dem fiir jedes 
Volumenelement 9 die thermischen Kopplungen mit seiner Umge- 

10 
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bung beriicksichtigt sind, Jedes Voluiaenelement 9 ist dabei - 
siehe FIG 3 - quaderformig ausgebildet. Es weist also drei 
Grundabmess-ungen A, B, C auf , die in der Kegel parallel zur 
Richtiing der Bandgeschwindigkeit v, zur Bandbreite b und zur 
5 Banddicke d orientiert sind. 

In einer ersten Ausgestaltung des erf indungsgemaSen Model- 
lierverf ahrens wird dem Rechner 4 gemaS FIG 4 in einem 
Schritt SI ein Erstzustand Z zugefuhrt. Der Erstzustand Z um- 
10 fasst fur jedes Volumenelement 9 des Stahlbandes 1 zunachst 

lokale Anteile pi, p2, p3 von modellierten Phasen des Stahls. 
Die Phasen konnen insbesondere Ferrit (Anteil pi) , Zementit 
(Anteil p2) und Austenit (Anteil p3 ) sein, 

15 Der Erstzustand Z umfasst weiterhin fiir jedes Volumenelement 
9 des Stahlbandes 1 eine GroSe H, die einen lokalen Energie- 
inhalt des Stahls des Volumenelements 9 beschreibt, Bei- 
spielsweise kann diese GroSe H die Enthalpie H des Volumen- 
elements 9 sein, Alternativ kommen auch die Tempera tur oder 

20 die Entropie in Frage. 

SchlieSlich umfasst der Erstzustand Z fiir mindestens eines 
der Volumenelemente 9 des Stahlbandes 1, vorzugsweise fiir je- 
des Volumenelement 9 des Stahlbandes 1, noch eine lokale Kon- 
zentrationsverteilung K mindestens eines im Stahl beweglichen 
Legierungselements . Das bewegliche Legierungselement kann 
insbesondere Kohlenstoff sein, Alternativ oder zusatzlich 
kommt z. B. auch Stickstoff in Frage. 

3 0 Dem Rechner 4 wird im Rahmen des Schrittes SI weiterhin zu- 
mindest eine gewiinschte EndgroSe f* zugefuhrt ♦ Gegebenen- 
falls konnen dem Rechner 4 im Rahmen des Schrittes SI auch 
noch gewunschte ZwischengroBen zugefuhrt werden, so dass dem 
Rechner 4 gegebenenf alls sogar ein gewunschter zeitlicher 

35 Verlauf der GroSe vorgegeben werden kann. 
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GemaS FIG 4 wird nun in einem Schritt S2 zunaclist ein Ein- 
flussgroSenverlauf festgelegt. Der Einf lussgroSenverlauf um- 
fasst fiir eine Vielzah.1 von zeitlich iinitiittelbar auf einander- 
folgenden Zeitpunkten jeweils eine Einf lussgroSe W. Die Ein- 
flussgroSe W entspricht beispielsweise der auf das Stahlvolu- 
men 1 (also das Stahlband 1 bzw. dessen betrachteten Bereich) 
aufzubringenden Kiihlmittelmenge unter Berucksichtigxing ande- 
rer Einfliisse auf das Stahlvolumen 1, wie beispielsweise von 
Transportrollen, der Warmekonvektion, der Warmeabstrahlung 
usw, . Sie umf asst fiir die Vielzahl von Flachenelementen 10 
des Stahlvolumens 1 (siehe FIG 2) jeweils einen lokalen Ein- 
fluss. Die lokalen Einflusse wirken dann iiber das jeweilige 
Flachenelement 10 auf das Stahlvolumen 1 ein. 

Gemafi FIG 4 wird nunmehr in einem Schritt S3 ein Anfangszu- 
stand ZA des Stahlvolumens 1 gleich dem Erstzustand Z und 
eine Zeitbasis t auf einen Anf angszeitpiinkt tO gesetzt. Der 
Anfangszeitpunkt tO entspricht dabei in der Kegel dem Zeit- 
punkt, zu dem das betrachtete Stahlvolumen 1, also z. B. ein 
Abschnitt des Stahlbandes 1, in die Kiihlstrecke 2 eintritt* 

Sodann wird in einem Schritt S4 anhand des Einf lussgroSenver- 
laufs die zu dem durch die Zeitbasis t bestimmten Zeitpunkt 
aufzubringende Einf lussgroSe W ermittelt. In einem Schritt S5 
eimittelt der Rechner 4 dann anhand des momentanen Anfangszu- 
stands ZA des Stahlvolumens 1 \xnd der iiber die Volumenober- 
f lache auf das Stahlvolumen 1 einwirkenden momentanen Ein- 
flussgroSe W einen Folgezustand ZF des Stahlvolumens 1. Der 
Rechner 4 lost dabei eine Warmeleitungsgleichung und eine 
Phaseniimwandlungsgleichung. Der Folgezustand ZF umf asst die- 
selben Elemente K, pi, p2, p3 , H wie der Anf angszustand ZA. 

Sodann wird in einem Schritt S6 die Zeitbasis t um eine zeit- 
liche Schrittweite At inkrementiert . Als nachstes wird in 
einem Schritt S7 gepruft, ob die Zeitbasis t einen Endzeit- 
punkt tl erreicht hat. Der Endzeitpunkt tl entspricht dabei 
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in der Kegel dam Zeitpunkt, zu dem das betrachtete Stahlvolu- 
men 1 die Kiihlstrecke 2 wieder verlasst. 

Wenn der Endzeitpunkt tl noch nicht erreicht ist, wird in ei- 
5 nem Schritt 38 der Anf angszustand ZA gleich dem iinmittelbar 
zuvor ermittelten Folgezustand ZF gesetzt und sodaiin zum 
Schritt S4 zuriickgesprungen . 

Wenn hingegen der Endzeitpunkt tl erreicht wurde, wird die 
aus den Schritten S4 bis S8 bestehende Schleife verlassen. 
Wenn danach keine weiteren Mafinahmen ergrif fen werden, wird 
nur ein sogenannter Prozessbeobachter realisiert. Vorzugs- 
weise aber wird, wie in FIG 4 dargestellt, in einem Schritt 
S9 anhand des nunmehr ermittelten Folgezustands ZA eine er~ 
wartete EndgroSe f , z. B. die Temperatur oder die Material- 
harte, ermittelt und mit der gewiinschten EndgroSe f vergli- 
chen. Wenn die erwartete EndgroSe f mit der gewiinschten End- 
groSe nicht korrespondiert/ wird - die entsprechende Re- 
chenleistung vorausgesetzt - in einem Schritt SIO der Verlauf 
der Einf lussgroSe W variiert. Sodann wird zum Schritt S3 zu- 
riickgesprungen. Anderenfalls ist der erf orderliche Kiihlmit- 
telmengenverlauf ermittelt, so dass der Rechner 4 nunmehr die 
Kiihlstrecke 2 in einem Schritt Sll entsprechend ansteuern 
kann, 

Der Vollstandigkeit halber sei dabei erwahnt, dass der Rech- 
ner 4, wenn ihm auSer der gewiinschten EndgroSe f* auch ge- 
wiinschte ZwischengroSen vorgegeben warden, selbstverstandlich 
den Kiihlmittelmengenverlauf erst bis zu der ersten gewiinsch- 
3 0 ten ZwischengroSe ermittelt, dann bis zur zweiten gewiinschten 
ZwischengroSe usw. , bis der gesamte Kiihlmittelmengenverlauf 
bis zur gewiinschten EndgroSe f* ermittelt ist, 

GemaS den FIG 1 und 4 wird das obenstehend beschriebene Mo- 
35 dellierverfahren online und in Echtzeit ausgefiihrt. Es ist 

aber selbstverstandlich ohne Weiteres auch eine Offline-Aus- 
fiihrung moglich. Dies ist in FIG 1 dadurch angedeutet, dass 
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die Verbindung der Kuhlstrecke 2 mit dem Rechner 4 als offen- 
bare Verbindung angedeutet ist, also nicht unmittelbar eine 
Ansteuerung der Kuhlstrecke 2 erfolgen muss. Weiterhin ist 
aus diesem Grund der Schritt Sll in FIG 4 nur gestrichelt 
5 eingezeichnet . 

Die Vorgehensweise gemaS FIG 4 erfordert einen sehr hohen Re- 
chenaufwand. Wenn die zur Ausfiihrung des Verfahrens gemaB FIG 
4 erf orderliche Rechenleistung nicht zur Verfiigung steht, 
10 aber dennoch eine Online- St euerung unter Nutzung des erfin- 
dungsgemaSen Modellierverf ahrens erfolgen soil, wird online 
und in Echtzeit ein Verfahren ausgefiihrt, das nachfolgend in 
Verbindung mit FIG 5 naher beschrieben wird. 

15 GemaS FIG 5 nimmt der Rechner 4 in einem Schritt S12 analog 
zum Schritt SI den Erstzustand Z entgegen. In einem Schritt 
S13 setzt der Rechner 4 analog zum Schritt S3 den Anfangszu- 
stand ZA gleich dem Erstzustand Z und die Zeitbasis t auf den 
Zeitpunkt tO • Sodann ermittelt der Rechner 4 in einem Schritt 

20 S14 eine gewiinschte FolgegroSe f* oder nimmt sie entgegen. 

In einem Schritt S15 ermittelt der Rechner 4 analog z-um 
Schritt S9 anhand des Anf angszustands ZA eine AnfangsgroSe f . 
Anhand der im Schritt S15 ermittelten AnfangsgroSe f und der 
gewiinschten FolgegroSe f * ermittelt der Rechner 4 dann in ei- 
nem Schritt S16 die Einf lussgroSe W. In einem Schritt S17, 
der mit dem Schritt Sll korrespondiert , steuert der Rechner 4 
schlieSlich die Kuhlstrecke 2 entsprechend der ermittelten 
Einf lussgroSe W an. Das Ansteuern erfolgt dabei, wie in FIG 5 
dargestellt, vorzugsweise sofort. Gegebenenf alls kann aber 
auch ein Ansteuern im nachsten Zyklus erfolgen. 

Die nachfolgenden Schritte S18 bis S21 korrespondieren mit 
den Schritten S5 bis S8 von FIG 4. Beziiglich dieser Schritte 
35 S18 bis S21 wird daher von einer detaillierten Erlauterung 
abgesehen. 
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Das Verfahren gemaS FIG 5 wird vorzugsweise z. B. bei einer 
Regelung eines die Kiihlstrecke 2 durchlauf enden Stahlbandes 1 
auf einen Verlauf der gewunschten FolgegroSe eingesetzt. 
In diesem Fall muss das obenstehend in Verbindung mit FIG 5 
beschriebene Verfahren selbstverstandlich individuell fiir je- 
den einzelnen Abschnitt des Stahlbandes 1 ausgefiilirt werden. 
Audi muss fiir jeden Abschnitt jeweils diejenige der Aufbrin- 
gungseinrichtungen 5 angesteuert werden, in deren Einflussbe- 
reich sich der betrachtete Abschnitt gerade befindet. Dies 
wird durch die eingangs erwahnte Wegverf olgung gewahrleistet . 
Das Verfahren gemaS FIG 5 wird in diesem Fall ferner parallel 
fiir alle in der Kiihlstrecke 2 befindlichen Abschnitte des 
Stahlbandes 1 ausgefiihrt. 

Die jeweils gewiinschte FolgegroSe f* kann dem Rechner 4 ex- 
plizit vorgegeben werden. Es ist aber auch moglich, dass dem 
Rechner 4 lediglich die gewiinschte FolgegroSe f* z. B. beim 
Austritt des Stahlbandes 1 aus der Kiihlstrecke 2 oder beim 
Erreichen eines Haspels, auf den das Stahlband 1 aufgehaspelt 
wird, vorgegeben wird. In diesem Fall ermittelt der Rechner 4 
selbsttatig anhand einer vorbestimmten Ermittlungsvorschrif t 
die gewiinschten FolgegroSen f^ fiir die einzelnen Iterationen, 

Die Regelung gemaB FIG 5 kann selbstverstandlich auch ander- 
weitig strukturiert werden, z. B. so, wie sie in der alteren 
DE 103 21 792.4 beschrieben ist. 

Wie erwahnt, werden in den Schritten S5 und S18 eine Warme- 
leitungsgleichung und eine Phasenumwandlungsgleichung gelost. 

Ein Beispiel einer Warmeleitungsgleichxang und einer Phasenum- 
wandlungsgleichung sind in FIG 6 dargestellt. Dort ist ein 
dreidimensionaler Ansatz fiir die Warmeleitung getroffen. Bei 
Bandern kann dieser Ansatz in aller Regel aber auf den eindi- 
mensionalen Ansatz gemaS FIG 7 vereinfacht werden, da der 
Warmef luss in Bandlangsrichtung und Bandquerrichtung in der 
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Regel vemachlassigbar ist. a/3x ist dabei die Ortsableitung 
in Banddickenrichtung. 

Ersichtlich. sind die Warmeleitungsgleicliung und die Phasenum- 
wandlungsgleichiing miteinander gekoppelt. Zur Losung der War- 
meleitungsgleichung und der Phasen-umwandlungsgleichung (= 
Schritt S5, S18) wird daher gemaS FIG 8 wie folgt vorgegan- 
gen: 

In einem Schritt S22 wird zunachst die lokale Temperatur T 
ermittelt. Die Ermittlung erfolgt anhand der Entlialpie H und 
der Anteile pi, p2, p3 der Phasen Ferrit, Zementit und Auste- 
nit. Auf Grxxnd der geringen GroSe des betrachteten VolTimen- 
eleitients 9 kann namlich angenoxnmen werden, dass die lokale 
Temperatur T innerhalb des Volumenelements 9 konstant ist. 
Damit weisen auch die Phasen des Stahls innerhalb dieses Vo- 
lumenelements 9 dieselbe lokale Temperatur T auf • Die Enthal- 
pie H des Volumenelements 9 lasst sich somit als 

H = p\m(T) + p2H2(T) + p3H3(T) ( 1 ) 

schreiben. Da ferner fiir jede der Phasen die Enthalpien Hi, 
H2 bzw. H3 der jeweiligen Phase als Funktion der lokalen Tem- 
peratur T - gegebenenfalls bei Austenit unter Berucksichti- 
gung des Kohlenstof f anteils - eineindeutig bestimmt sind, 
lasst sich durch die Gleichung 1 die lokale Temperatur T ohne 
Wei teres ermi t teln . 

Fiir die im Schritt S22 ermittelte lokale Temperatur T werden 
dann in einem Schritt S23 - getrennt fiir die Phasen Ferrit, 
Zementit und Austenit - die Gibbsschen freien Enthalpien Gl 
bis G3 der einzelnen Phasen als Funktion des Anteils des be- 
weglichen Legierungselements ermittelt. Beispiele derartiger 
Verlaufe sind in den FIG 9 xmd 10 dargestellt- Bei dem Bei- 
spiel gemaS FIG 9 liegt die lokale Temperatur T dabei ober- 
halb der Umwandlizngs temperatur von Austenit in Per lit, bei 
dem Beispiel gemaS FIG 10 darunter. 
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Es ist moglich/ stets anzunehmen , dass alle drei Phasen Fer- 
rit, Zementit und Austenit vorliegen und erst spater - siehe 
die nachfolgenden Erlauterungen zu FIG 18 - zu entscheiden, 
welche Phasen vorliegen und welche PhaseniamwandlTongen erfol- 
gen. GemaS FIG 8 wird aber zxinachst in einem Schritt S24 die 
Anzahl der zu betrachtenden Phasen ermittelt. Dies geschieht 
wie folgt: 

Zunachst wird das System Ferrit-Austenit heraus gegriffen. 
Fiir dieses System wird iiberpriift, ob es bei dem im Vol\imen- 
element 9 vorhandenen Gesamtanteil des beweglichen Legie- 
rungselements rein ferritisch, rein austenitisch oder ge- 
mischt stabil ist und welche Phasenverteilung gegebenenf alls 
im stabilen Zustand vorliegt. Der Gesamtanteil des bewegli- 
chen Legierungselements kann anhand der Konzentrationsvertei- 
lung ohne Weiteres ermittelt werden. 

Die Ermittlujig der stabilen Phase (n) und gegebenenf alls der 
Phasenverteiluxig geschieht dadurch, dass versucht wird, ein 
Minimum fiir die gesamte Gibbssche freie Enthalpie G eines 
derartigen Systems zu ermitteln. Die konkrete Vorgehensweise 
ergibt sich dabei so, wie sie in den Seiten 16 bis 18 der al- 
teren deutschen Patentanmeldung 102 51 716.9 beschrieben ist. 
Diese Vorgehensweise liefert zugleich auch Konzentrationen kl 
xind k3, in denen das bewegliche Legierungselement - typi- 
scherweise Kohlenstoff - an einer etwaigen Phasengrenze 11 
zwischen Ferrit und Austenit im ferritischen bzw. im austeni- 
tischen Bereich gegebenenf alls vorliegt - 

In analoger Weise werden sodann die Mischungssysteme Ferrit- 
Zementit (= Perlit) und Austenit-Zementit vmtersucht. Die Un- 
tersuchung des Mischungssystems Austenit-Zementit liefert da- 
bei zugleich auch Konzentrationen k2 und k4, in denen das be- 
wegliche Legierungselement an einer etwaigen Phasengrenze 12 
zwischen Zementit und Austenit im zementitischen bzw. im 
austenitischen Bereich gegebenenf alls vorliegt. 
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Die drei fiir die drei Zweiphasen-Systeme ermittelten minima- 
len Gibbsschen freien Enthalpien G werden miteinander vergli- 
chen iind das Zweiphasen-System mit der insgesamt minimalen 
Gibbsschen freien Enthalpie G herausgegrif f en. Wenn dieses 
Zweiphasen-System Austenit enthalt, lasst sich anhand der 
Gibbsschen freien Enthalpien ermitteln, ob der stabile Zu- 
stand eine oder zwei Phasen umfasst, welche Phase dies gege- 
benenfalls ist bzw. welche Phasen dies gegebenenf alls sind 
und welche Konzentrationen kl, k3 bzw. k2, k4 des beweglichen 
Legieriingselements beidseits der Phasengrenze 11 von Austenit 
zu Ferrit bzw. der Phasengrenze 12 von Austenit zu Zementit 
jeweils vorliegen. 

Aus dem Anf angszustand ZA ist femer bekannt, welche Phasen 
anfanglich zu welchen Anteilen pi bis p3 vorliegen. 

Wenn das Zweiphasen-System mit deir minimalen Gibbsschen 
freien Enthalpie G Austenit enthalt, lasst sich also durch 
eine Zusammenschau der anfanglich vorliegenden Phasen und der 
als stabil ermittelten Phasen ermitteln, ob ein Reinsystem 
mit nur einer der Phasen Ferrit, Austenit uxid Zementit oder 
ein Mischsystem mit zwei oder sogar drei der Phasen Ferrit, 
Austenit und Zementit betrachtet werden muss. Auf Grund des 
Umstands, dass anhand der Gibbsschen freien Enthalpien G, Gl 
bis G3 auch ermittelbar ist, welche Phase stabil ist bzw. 
welche Phasen stabil sind, ist dariiber hinaus auch die Rich- 
tiing einer etwaigen Phasenumwandlung bereits bekannt. 

Wenn das Zweiphasen-System mit der insgesamit minimalen 
Gibbsschen freien Enthalpie G das Perlitsystem ist und der 
Anfangszustand ZA kein Austenit enthalt (p3 = 0), liegt ein 
vollstandig umgewandeltes Gefiige, also ein Zweiphasen-System 
Ferrit-Zementit vor, das keiner Phasen^lmwandlung mehr vmter- 
liegt. Wenn hingegen der Anfangszustand ZA Austenit enthalt 
(p3 > 0) und das Perlitsystem die minimale Gibbssche freie 
Enthalpie G aufweist, lasst sich nicht ohne Weiteres ent- 
scheiden, ob nun eine oder beide der Phasen Ferrit und Zemen- 
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tit gebildet warden. In diesem Fall wird daher vorlaufig an- 
genommen, dass beide Phasen gebildet werden, also ein Drei- 
phasensystem Ferrit-Zementit-Austenit zu betrachten ist. 

In einem Schritt S25 wird als nSchstes gepriift, ob in dem be- 
trachteten Volinaenelement 9 mehr als eine Phase vorliegt. 
Wenn dies nicht der Fall ist, erfolgt selbstverstandlich 
keine Phasenumwandlung. Dennoch wird in einem Schritt S26 fiir 
dieses Volvunenelement 9 eine Dif fusionsgleichiang fiir das be- 
wegliche Legieriingselement angesetzt wcid gelost, um etwaige 
Konzentrationsschwankungen des beweglichen Legierungselements 
innerhalb des betrachteten Volumene laments 9 auszugleichen. 
Diese Vorgehensweise ist insbesondere bei Austenit, bei dem 
der Gahalt an Kohlenstoff stark schwanken kann, von Bedeu- 
tiing. Der Ansatz und die Losiong einer derartigen Diffusions- 
gleichung ist Fachleuten allgemein bekannt. Auf den Ansatz 
der Diffusionsgleichung und deren Losung wird daher nachste- 
hend nicht naher eingegangen. 

Wenn in dem betrachteten Voliamenelement 9 mehr als eine Phase 
vorliegt, wird als nachstas in einem Schritt S27 abgepriift, 
ob alle drei Phasen vorliegen. Wenn dies nicht der Fall ist, 
liegt ein Zweiphasen-System vor. In diesem Fall wird in einem 
Schritt S28 uberpriift, ob eine der beiden vorliegenden Phasen 
Austenit ist. Wenn dies nicht der Fall ist, liegt ein voll- 
standiges Perlit-System vor, also eine Schichtstruktur aus 
Ferrit und Zementit- Diese Struktur ist im Wesentlichen sta- 
bil. In diesem Fall mussen daher kaina weiteren MaSnahmen er- 
griffen werden. 

Wenn hingegen eine der Phasen Austenit ist und somit die an- 
dere Phase Ferrit oder Zementit, existiert eine Phasengrenze 
11 zwischen Austenit und Ferrit (siehe FIG 11 und 12) oder 
eine Phasengrenze 12 zwischen Austenit und Zementit (siehe 
FIG 11 und 13 ) . In diesem Fall wird als vereinf achende MaS- 
nahme angenommen, dass die Phasengrenze 11 bzw. 12 als Recht- 
eck ausgebildet ist, das eine Langssaite und eine Querseite 
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aufweist. Die Langsseite und die Querseite des Rechtecks kor- 
respondieren dabei mit den Grundabmessungen A und B des be- 
trachteten Volumenelements 9. Verschiebungen der Phasengrenze 
11, 12 erfolgen in diesem Fall, wie insbesondere aus FIG 11 
ersichtlich ist, parallel zur Grundabmessung C. 

Wie Fachleuten allgemein bekannt ist und auch aus FIG 12 er- 
sichtlich ist, ist in Ferrit nur sehr wenig Kohlenstoff ent- 
halten. An der Phasengrenze 11 zum austenitischen Bereich 
hingegen erfolgt ein Konzentrationssprung. Die Konzentratio- 
nen kl, k3 , in denen das bewegliche Legierungselement an der 
Phasengrenze 11 vorliegt, sind dabei die Konzentrationen kl, 
k3, die obenstehend im Schritt S24 fur das Ferrit-Austenit- 
System ermittelt wurden, 

Wenn, wie in FIG 12 dargestellt, eine Umwandlung von Austenit 
in Ferrit erfolgen soil, muss der ^^Konzentrationsberg" des 
beweglichen Legierungselements im austenitischen Bereich in 
der Nahe der Phasengrenze 11 abdif fundieren. Umgekehrt muss, 
wenn eine Umwandlung von Ferrit in Austenit erfolgen soil, 
das ^^Konzentrationstal'' des beweglichen Legierungselements im 
austenitischen Bereich in der Nahe der Phasengrenze 11 stan- 
dig aufgefullt werden. Die numerische oder analytische Losixng 
dieser Aufgabe ist allgemein als Stefan-Problem bekannt. Sie 
erfolgt in einem Schritt S29. Das Stefan-Problem lasst sich 
im vorliegenden Fall wie folgt formulieren: 

t XX 

D ist dabei die - gegebenenf alls temperaturabhangige - Diffu- 
sionskonstante des beweglichen Legierungselements im Auste- 
nit. K' ist die Konzentration K. Die Indizes t und x bedeuten 
die Ableitung nach der Zeit bzw. dem Ort in Verschieberich- 
tung der Phasengrenze 11 . Die. Verschieberichtung der Phasen- 
grenze 11 muss dabei nicht notwendigerweise die Banddicken- 
richtung sein. 
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Fiir die Gleichung 2 sind die Aaf angsbedingungen durch die 
ortliche Konzentrationsverteilung K' des beweglichen Legie- 
riingseleitients im Anf angszustand ZA definiert, Beim Losen des 
Stefan-Problems mussen dann nur noch die Randbedingungen be- 
trachtet werden, dass im f erristischen bzw. im austenitischen 
Bereich an der Phasengrenze 11 die Konzentrationen kl bzw. k3 
vorliegen, dass das bewegliche Legierungselement das betref- 
fende Volumenelement 9 nicbt verlassen kann und dass sich die 
Verschiebung 6x' der Phasengrenze 11 gemafi der Stef an-Bedin- 
gung zu 



ergibt, wobei 6t die beim Losen des Stefan-Problems benutzte 
zeitliche Schrittweite ist, Diese kann dabei ein Bruchteil 
(1/2, 1/3, 1/4 •..) der zeitlichen Schrittweite At oder aber 
gleich der zeitlichen Schrittweite At sein, 

Durch Losen des Stef an- Problems im Schritt S29 wird also er- 
mittelt, ob und wie sich die Konzentrationsverteilung K bzw. 
K' des beweglichen Legierungselements im austenitischen Be- 
reich des betrachteten Voliimenelements 9 andert. Zugleich 
wird dabei auch ermittelt, ob und um welches AusmaS 5x' sich 
die Phasengrenze 11 dadurch verschiebt. 

Falls die weitere Phase nicht Ferrit, sondern Zementit ist, 
ergibt sich ein prinzipiell gleich gearteter Losungsansatz 
und auch eine prinzipiell gleich geartete Losung. Der Unter- 
schied ist nur, dass bei der Bildung von Zementit im Zementit 
Kohlenstoff stark angereichert ist und im austenitischen Be- 
reich in der Nahe der Phasengrenze 12 zum Zementit Kohlen- 
stoff verarmt ist. Auch dies ist Fachleuten allgemein bekannt 
und in FIG 13 dargestellt. Es andern sich also die Randbedin- 
gungen. Die VerschiebTing 5x" der Phasengrenze 11 gemafi der 
Stefan-Bedingung ergibt sich in diesem Fall zu 
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K" ist dabei wieder die Konzentration K. 

Wenn hingegen im Schritt S28 festgestellt wurde, dass alle 
drei Phasen vorliegen bzw. vorliegen konnten, ergibt sich ein 
komplizier teres Problem. Derm dann muss der Aasatz gewahlt 
werden, dass eine Schiclitstruktur existiert, die - siehe FIG 
14 und 15 - abwechselnd aus einer Ferritschicht 13 und einer 
Zementitschicht 14 besteht. Die Schichtstruktur grenzt an ei- 
nen austenitischen Bereich. 15 an. Die FIG 16 und 17 zeigen 
jeweils einen Verlauf der Konzentration K' bzw. K" des beweg- 
lichen Legierungselements in einer Ferritschicht 13 bzw. ei- 
ner Zementitschicht 14 und im austenitischen Bereich 15 vor 
diesen Schichten 13 , 14 . 

Im Rahmen des Modells wird dabei nach wie vor angenommen, 
dass die Phasengrenzen 11, 12 als Rechtecke ausgebildet sind, 
die jeweils eine Langsseite und eine Querseite aufweisen. Die 
Langsseiten korrespondieren weiterhin mit der ersten Grundab- 
messung A. Die Querseiten verlauf en nach wie vor parallel zur 
zweiten Grundabmessung B. Auch wird nach wie vor angenommen, 
dass Verschiebungen der Phasengrenzen 11, 12 parallel zur 
dritten Grundabmessung C erfolgen. 

Fiir jede der Phasengrenzen 11 und 12 wird nunmehr in einem 
Schritt S30 ein Stefan-Problem angesetzt und gelost. Der 
Schritt S3 0 ist in FIG 18 detaillierter dargestellt. 

Das Stefan-Problem fiir die Phasengrenze 11 zwischen Ferrit 
und Austenit gehorcht gemaS FIG 18, siehe dort Schritt S31, 
folgender GesetzmaSigkeit : 

K' ist dabei die Konzentration des beweglichen Legierungsele- 
ments vor der Phasengrenze 11. 
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Das Stefan- Problem fur die Phasengrenze 12 zwischen Zementit 
und Austenit gehorcht folgender GesetzmaSigkeit : 



(6) 



10 




K" ist dabei - analog zu K' - die Konzentration des bewegli- 
clien Legieriingselements vor der Phasengrenze 12. 

LI und L2 sind Kopplungsterme. Sie sind - siehe FIG 14 - 
Fxonktionen eines Lamellenabstands 1 der Schichtstruktur und 
eines Antells q der Zementitphase an der Schichtstruktur so- 
wie dessen zeitlicher Ableitung. Ausgeschrieben konnen sie z 
B . lauten : 
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Ll = 



1 



/2 t 



xond 



(7) 



L2 = 



1 ( 48(1 - q)D ..'^ 



(8) 




Wie aus den Gleichungen 5 bis 8 ersichtlich ist, sind die 
20 beiden Stef an-Probleme auf Grund der Kopplungsterme Ll, L2 

auf der rechten Seite der Gleichungen 5 und 6 miteinander ge- 
koppelt. 

Auch in den Gleichungen 5 bis 8 stehen die Indizes t und x 
25 wieder fur die Ableitiing nach der Zeit bzw. dem Ort in Ver- 
schieberichtting der Phasengrenze 11 bzw. 12. 

Fiir die Verschiebungen 5x' , 5x" der Phasengrenzen 11, 12 gel- 
ten weiterhin die Stef an-Bedingungen (siehe Gleichungen 3 und 
30 4). Es ist daher moglich, die Verschiebungen 5x' , 8x" gemaB 
den Gleichungen 3 und 4 in einem Schritt S3 2 gleichzusetzen. 
Dann ist der Anteil q derart zu bestimmen, dass die Verschie- 
bung 5x' der Phasengrenze 11 gemaS Gleichixng 3 und die Ver- 
schiebung 5x" der Phasengrenze 12 gemaS Gleichung 4 denselben 
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10 





Wert annehjnen. Der Anteil q ist in diesem Fall also derart zu 
bestimmen, dass die Phasengrenzen 11, 12 stets nebeneinander 
angeordnet bleiben. 

Die Gleichlieit der Verschiebungen 5x' , 5x" der Phasengrenzen 
11, 12 kann also dadurch erreicht werden, dass der Anteil q 
der Zementitphase an der Schichtstruktur des Perlit geeignet 
gewahlt wird. Man kann also - z. B. durch. Einsetzen imd Aus- 
rechnen (trial and error) - ermitteln, bei welcliem Anteil q 
der Zementitphase die Verschiebungen 5x' , 5x" der Phasengren- 
zen 11, 12 iibereinstimmen. 

« 

Anliand des so ermittelten Anteils q lasst sich. dann auch 
bestimmen, ob tatsachlich Perlit gebildet wird oder ob nur 
15 eine der beiden Phasen Ferrit und Zementit gebildet wird. 

Denn nur wenn der Anteil q zwischen Null und Eins liegt, wird 
Perlit gebildet. Ist der Anteil q hingegen groSer als Eins, 
wird ausschlieSlich Zementit gebildet. Ist er umgekehrt klei- 
ner als Null, wird ausschlieSlich Ferrit gebildet. 

20 

Der Anteil q kann bei dieser Vorgehensweise aus numerischen 
Griinden beliebige Werte annehmen, also insbesondere auch 
kleiner als Null oder groSer als Eins sein. Diese Werte sind 
aber physikalisch sinnlos. Der Anteil q wird daher in Sclirit- 
ten S3 3 bis S3 6 gegebenenf alls entsprecbend korrigiert. 

Die Vorgehensweise gemaS Schritt 32 kann zu niimerischen Prob- 
lemen fuhren. Alternativ ist es daher moglich, unterschiedli- 
che Verschiebungen 8x' , 5x" zuzulassen. In diesem Fall wird 
3 0 der Anteil q des Zementit an der Schichtstruktur gemaS 

"''"-la- ' ' 

in Schritten S32a land S32b derart nachgefiihrt, dass die Lagen 
35 der Phasengrenzen 11, 12 aufeinander zustreben. Die Schritte 
S32a und S32b werden dabei gemaE FIG 18 anstelle des Schrit- 
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tes S3 2 ausgefiihrt. a ist eine geeignete Proper tionalitats- 
konstante. Ihr Wert ist grofier als Null. 

In diesem Fall wird eine mittlere Verschiebung 5x der Phasen- 
grenzen 11, 12 zu 

def iniert • 

Unabhangig davon, welche dieser beiden Vorgehensweisen 
(Scliritt S32 oder Schritte S32a, S32b) ergriffen wird, wird 
aber in beiden Fallen erst durch die Bestimmung des Anteils q 
ermittelt und abgeleitet, welche Phasen in welchem Umfang ge- 
bildet werden. 

Wenn Per lit gebildet wird, der Anteil q also zwischen Null 
und Eins liegt, muss noch der Lamellenabstand 1 bestimmt wer- 
den. Dies geschieht gemaS FIG 19 wie folgt: 

Mit der Verschiebung 5x der Phasengrenzen 11, 12 und dem An- 
teil q ist bekannt, welche Phasenumwandlungen erfolgen und in 
welchen Anteilen sie erfolgen. Die Phasenumwandlungen liefern 
einen Beitrag 5e1 zur Energiebilanz . Der Beitrag 5e1 ist von 
dem Vo lumen abhangig, in dem die Phasenumwandlungen erfolgen. 
Es gilt also 

P ist dabei eine vorab ermittelbare Proportionalitatskon- 
stante. 

Es ergibt sich femer eine Anderung der Flache einer Grenz- 
schicht 16 zwischen Ferrit und Zementit. Auch die Anderung 
dieser Flache liefert einen Beitrag 5E2 zur Energiebilanz. 
Dieser Beitrag 5e2 ist proportional zur Anderung der Flache 
der Grenzschicht 16. Es gilt also 
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m2 = 2yAdx (12) 

7 ist dabei wieder eine vorab ermittelbare Proportionalitats- 
konstante. Der Faktor 2 ergibt sich dadurch, dass pro Phasen- 
5 grenze 11, 12 bzw. pro Schicht 13, 14 je eine Grenzschicht 16 
vorhanden ist und zwei Phasengrenzen 11, 12 bzw. zwei Schich- 
ten 13, 14 betrachtet werden. 



10 




15 



20 




30 



Anhand der Gleichxangen 11 und 12 lasst sich somit ein kriti- 
scher Lamellenabstand 1' ermitteln, bei dam die Energiebi- 
lanz, welche die beiden Beitrage 5e1 tind 5e2 berucksichtigt, 
ausgeglichen ist. Fiir diesen kritischen Lamellenabstand I' 
gilt 

Der Lamellenabstand I, also die S\unme der Schichtdicken 11 
und 12, wird nun auf das ca. eineiiihalb- bis dreifache dieses 
kritischen Lamellenabstand 1' gesetzt, z. B. auf das ca. 
zweifache. Die Schichtdicken 11 bzw. 12 ergeben sich dann zu 



/I = (1 - q)l Tond 



(14) 



I2 = ql (15) 

Die Schichtdicken 11, 12 sind proportional zu FlachenmaSen 
Fl, F2, welche die Phasengrenzen 11, 12 aufweisen. Je eines 
der FlachenmaSe Fl, F2 ist in FIG 14 durch Schraffieren ent' 
spr echend hervorgehoben . 



Wenn in einem der Schritte S29, S30 eine Verschiebung 6x' , 
5x" einer Phasengrenze 11, 12, ermittelt wurde, hat sich so- 
mit die Lage der Phasengrenze 11, 12 geandert. In diesem Fall 
wird in einem Schritt S37 (selbstverstandlich xanter Beachtung 
35 des Vorzeichens der mittleren Verschiebxing 8x) eine neue Lage 
der Phasengrenzen 11, 12 oder 11 und 12 bestimmt. In einem 
Schritt S3 8 wird dann anhand einer nicht linearen Fxinktion 
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der Lage der Phasengrenze 11, 12 bzw. der Phasengrenzen 11, 
12 der Anteil p3 des Austenits ermittelt. Die nicht lineare 
Funktion beriicksiclitigt dabei insbesondere den Umstand, dass 
das Stefan- Problem des Schrittes S29 bzw. die Stef an-Probleme 
5 des Schrittes S3 0 eindimensional angesetzt und gelost warden, 
in der Real i tat aber eine dreidimensionale Umwandlung er- 
folgt . 

In einem Schritt S3 9 werden sodann die Anderungen der Anteile 
10 pi, p2 der beiden anderen Phasen Ferrit und Zementit be- 
^ stimmt. Falls dabei die Schritte S37 bis S39 vom Schritt S29 
B aus erreicht werden, ist selbstverstandlich der Anteil q Null 

Oder Eins . 



15 Mit der mittleren Verschiebung 5x der Phasengrenzen 11, 12 

und dem Anteil q des Zementit an der Schichtstruktur (0 < q < 
1) lasst sich also auf einfache Weise ermitteln, welche Ande- 
rungen sich fiir die Anteile pi, p2 , p3 der Phasen des Stahls 
beziiglich des betrachteten Volumenelements 9 ergeben. 

20 

Trotz der vereinf achenden Annahmen des obigen Modellierver- 
fahrens ergibt sich ein erheblicher Rechenauf wand . Daher wer- 
den gemaS FIG 8 vor dem Schritt S22 Schritte S40 und S41 und 
^ nach dem Schritt S39 ein Schritt S42 eingefugt. 

Im Schritt S40 werden die Volumenelemente 9 zu Gruppen zusam- 
men gefasst. Beispielsweise kann eine Zusammenf assung mehre- 
rer Volumenelemente 9 erfolgen, die in Richtung der Band- 
breite b, der Banddicke d und/oder der Bandgeschwindigkeit v 

3 0 aneinander angrenzen, Auch Kombinationen hiervon sind mog- 

lich. Im Schritt S41 wird dann pro Gruppe von Volumenelemen- 
ten 9 je ein einziges der Volumenelemente 9 ausgewahlt, Nur 
fiir die ausgewahlten der Volumenelemente 9 werden in den 
Schritten S22 bis S3 9 die Dif f erentialgleichung, das Stefan- 

35 Problem bzw. die Stef an-Probleme gelost und die Anteile pi 
p2, p3 der Phasen rechnerisch ermittelt. 
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GemaS FIG 8 wird bei jeder Gruppe das Losungsergebnis des be- 
trachteten Voliimenelements 9 im Schritt S42 fur die anderen 
Voliomenelemente 9 der jeweiligen Gruppe ubernommen. Dies 
stellt die einfachste Vorgehensweise dar. Es ware aber auch 
denkbar, die Anteile pi, p2, p3 der Phasen von Volumenelemen- 
ten 9, deren Phasenverteilung nicht explizit berechnet worden 
ist, durch linear e oder nicht lineare Interpolation zu ermit- 
teln. 

Auf Grund dieses Umstands, namlich. dass bei jeder Gruppe von 
Volumenelementen 9 nur die Konzentrationsverteilung K eines 
der Volumenelemente 9 errechnet wird und fiir das Losen der 
Warmeleitungsgleichung nicht die Konzentrationsverteilung K, 
sondern nur die Anteile pi, p2 , p3 der vorhandenen Phasen be- 
notigt werden, kann gegebenenf alls auch die Konzentrations- 
verteilung K nur fiir dieses eine Volumenelement 9 jeder 
Gruppe vorgegeben werden. 

Die Warmeleitungsgleichung hingegen wird in einem Schritt S43 
individuell fiir jedes der Volumenelemente 9 gelost. Der 
Schritt S43 kaim dabei alternativ vor oder nach dem Losen der 
Phasenumwandlujigsgleichung in den Schritten S22 bis S42 aus- 
gefiihrt werden. 

« 

Das Losen der Warmeleitxingsgleichung ist ohne Weiteres mog- 
lich. Denn die Temperaturen T der einzelnen Voliiitienelemente 9 
sind - siehe obige Erlauterungen zum Schritt S22 - ohne Wei- 
teres eriaittelbar, so dass auch deren Ortsgradient sofort er- 
mittelbar ist. Da auch die Anteile pi, p2, p3 der Phasen der 
einzelnen Voliamenelemente 9 ebenfalls bekannt sind, lasst 
sich auch die Warmeleitf ahigkeit X der einzelnen Volumenele- 
mente 9 ohne Weiteres ermitteln. Da schlieSlich die Dichte p 
im Wesentlichen eine Konstante ist tind die Enthalpie H der 
einzelnen Volumenelemente 9 direkt gegeben ist, ist somit 
auch die Warmeleitungsgleichung als Ganzes losbar . 
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Mittels der vorliegenden Erf indung wird also ein physikali- 
schen Modell erstellt, bei dem auf Basis der Gibbsschen 
freien Enthalpie G und des Dif fusionsgesetzes (Stefan- Prob- 
lem) das Umwandlungsgescliehen hinsichtlich entstehender Pha- 
5 senanteile pi, p2 , p3 und Umwandlungsgeschwindigkeit mit sehr 
hoher Genauigkeit beschreibbar ist. Es werden allgemein gul- 
tige Ergebnisse ermoglicht, welche auch die Behandlung noch 
unbekannter Mater ialien und Werkstoffe erlauben. Die Ansatze 
sind daruber hinaus nicht nur zur Temperaturberechnung, son- 

10 dern auch zur Berechniing von Gefiigestruktur und KorngroBe an- 

^ wendbar • 

^ Das erf indungsgemaSe Modellierverf ahren kann an beliebiger 
Stelle eingesetzt werden, also insbesondere beispielsweise 

15 auch bei Kuhlvorgangen zwischen Walzgerusten einer WalzstraSe 
Oder bei sogenanntem ferritischem Walzen. Es eignet sich auch 
zur Beschreibung der Ruckumwandlung beim Aufheizen des Stah- 
les • 
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Patentanspruche 

1. Rechnergestiitztes Model lierverfahr en fiir das Verhalten ei- 
nes StahlvolTomens (1) xtiit einer Voliimenoberf lache, 

- wobei ein Rechner (4) anhand eines momentanen Anfangszu- 
stands (ZA) des Stahlvoliamens (1) und mindestens einer uber 
die volumenoberflache auf das Stahlvoliamen (1) einwirkenden 
momentanen Einf lussgroSe (W) durch Losen einer Warmelei- 
tTingsgleichung und einer Phaseniomwandlungsgleichung einen 
Folgezustand (ZF) des Stahlvoliamens (1) ermittelt, 

- wobei die mindestens eine Einf lussgroSe (W) fiir eine Anzahl 
von Flachenelementen (10) der Volumenoberflache jeweils ei- 
nen lokalen Einfluss umfasst und die lokalen Einfliisse iiber 
das jeweilige Flachenelement (10) auf das Stalilvolumen (1) 
einwirken, 

- wobei der Anf angszustand (ZA) und der Folgezustand (ZF) fiir 
eine Anzahl von Voliomenelementen (9) des Stahlvoliamens (1) 
jeweils lokale Anteile (pi, p2, p3) von modellierten Phasen 
des Stahls und eine einen lokalen Energieinhalt des Stahls 
beschreibende GroSe (H) umfassen, 

- wobei die modellierten Phasen des Stahls Austenit und eine 
erste weitere Phase umfassen, in die Austenit umwandelbar 
ist \xn.d die in Austenit umwandelbar ist, 

dadurch gekennzeichnet, 

- dass der Anf angszustand (ZA) und der Folgezustand (ZF) fiir 
mindestens eines der Volumenelemente (9) auch eine lokale 
Konzentrationsverteilung (K) mindestens eines im Stahl be- 
weglichen Legierungselements \imfassen, 

- dass im Rahmen der Umwandlxingsgleichung fiir das mindestens 
eine Volumenelement (9) ermittelt wird, welche Konzentrati- 
onen (kl, k3 ; k2 , k4) des mindestens einen beweglichen Le- 
giervmgselements beidseits einer ersten Phasengrenze (11, 
12) zwischen Austenit und der ersten weiteren Phase vorlie- 

gen, 

- dass durch Losen eines ersten Stefan-Problems ermittelt 
wird, ob und wie sich die KonzentrationsverteilTong (K) des 
mindestens einen beweglichen Legierungselements im austeni- 
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tischen Bereich des betrachteten Volumene laments (9) andert 
und ob und ixm welches AusmaS (5x, 5x' , 8x" ) sich. die erste 
Phasengrenze (11, 12) dadurch verschiebt, und 

- dass die lokalen Anteile (pi, p2, p3 ) der Phasen anhand ei- 
ner durch das AusitiaS (5x) der Verschiebung der ersten Pha- 
sengrenze (11, 12) bestiinmten Lage der ersten Phasengrenze 
(11, 12) ermittelt werden. 

2. Modellierverfahren nach Anspruch 1, 
dadurch gekennzeichnet, 

- dass die modellierten Phasen des Stahls auch eine zweite 
weitere Phase umfassen, in die Austenit uitiwandelbar ist und 
die in Austenit umwandelbar ist, 

- dass fiir das betrachtete Volumene lament (9) im Rahmen der 
Umwandlungsgleichung auch ermittelt wird, welche Konzentra- 
tionen (k2, k4; kl, k3 ) des mindestens einen beweglichen 
Legierungselements beidseits einer zweiten Phasengrenze 
(12, 11) zwischen Austenit und der zweiten weiteren Phase 
vorliegen, 

- dass durch zusatzliches Losen eines zweiten Stefan-Problems 
ermittelt wird, ob und wie sich die Konzentrationsvertei- 
lung (K) des mindestens einen beweglichen Legierungsele- 
ments im austenitischen Bereich des betrachteten Volumen- 
elements (9) andert und ob Tind um welches AusmaS (5x", 5x' ) 
sich die zweite Phasengrenze (12, 11) dadurch verschiebt, 

- dass die Stef an-Probleme miteinander gekoppelt sind, 

- dass den Phasengrenzen (11, 12) FlachenmaSe (Fl, F2) zuge- 
ordnet sind, 

- dass ein Anteil (q) des der zweiten Phasengrenze (12) zuge- 
ordneten FlachenmaSes (F2) an der Summe der FlachenmaSe 
(Fl, F2) bestimmt wird und 

- dass die lokalen Anteile (pi, p2, p3 ) auch vom Anteil (q) 
des der zweiten Phasengrenze (12) zugeordneten FlachenmaiSes 
(F2) an der Summe der FlachenmaSe (Fl, F2) abhangen. 

3. Modellierverfahren nach Anspruch 2, 

dadurch gekennzeichnet, dass der An- 
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tail (q) des der zweiten Phasengrenze (12) zugeordneten Fla- 
chenmaSes (F2) an der Siomme der FlachenmaSe (Fl, F2) derart 
bestimmt wird, dass die Phasengrenzen (11, 12) stets neben- 
einander angeordnet bleiben. 

4. Modellierverfahren nach Anspruch 2, 

dadurch. gekennzeichnet, dass der An- 
teil (q) des der zweiten Phasengrenze (12) zugeordneten Fla- 
chenmaSes (F2) an der Sunune der FlachenmaSe (Fl, F2) derart 
nachgefiilirt wird, dass die Phasengrenzen (11, 12) aufeinander 
zu streben. 

5. Model lierverfalir en nach Anspruch 2, 3, oder 4, 
dadurch gekennzeichnet, dass anhand 
des Anteils (q) des der zweiten Phasengrenze (12) zugeordne- 
ten FlachenmaSes (F2) an der Summe der FlachenmaSe (Fl, F2) 
abgeleitet wird, ob Austenit nur in die erste weitere Phase, 
nur in die zweite weitere Phase oder sowohl in die erste als 
auch in die zweite weitere Phase umgewandelt wird. 

6. Modellierverfahren nach einem der obigen Anspriiche, 
dadurch gekennzeichnet, 

- dass das betrachtete Volumenelement (9) quaderf ormig ausge- 
bildet ist und drei Grundabmessungen (A, B, C) aufweist, 

- dass die erste Phasengrenze (11, 12) als Rechteck mit einer 
ersten Langsseite und einer ersten Querseite ausgebildet 
ist und 

- dass die erste Langsseite mit einer ersten der Grundabmes- 
sungen (A, B, C) korrespondiert, die erste Querseite paral- 
lel zu einer zweiten der Grundabmessungen (A, B, C) ver- 
lauft und Verschiebungen (Sx', 8x" ) der ersten Phasengrenze 
(11, 12) parallel zu der dritten der Grundabmessungen (A, 
B, C) erfolgen. 

7. Modellierverfahren nach Anspruch 6 und einem der Anspriiche 
2 bis 5/ dadurch gekennzeichnet. 
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- dass die zweite Phasengrenze (12, 11) als Rechteck mit ei- 
ner zweiten Langsseite und einer zweiten Querseite ausge- 
bildet ist vind 

- dass die zweite Langsseite mit der ersten der Grundabmes- 
sungen (A, B, C) korrespondiert , die zweite Querseite pa- 
rallel zur zweiten der Griindabmessungen (A, B, C) verlauft 
und Verschiebungen (5x", 5x' ) der zweiten Phasengrenze (12, 
11) parallel zur dritten der Grundabmessungen (A, B, C) er- 
f olgen. 

8. Modellierverfahren nach Anspruch 7, 

dadurch gekennzeichnet, dass die 
Summe (1) der Querseiten der Phasengrenzen (11, 12) ungefahr 
gleich dem 1,5- bis 3-fachen eines kritischen Lamellenabstan- 
des (1') ist, bei dem eine Energiebilanz , die einerseits mit 
dem Verschieben der Phasengrenzen (11, 12) korrespondierende 
Phasenumwandlungen des Stahls und andererseits mit dem Ver- 
schieben der Phasengrenzen (11, 12) korrespondierende Ande- 
rungen der Flache einer Grenzschicht (16) zwischen der ersten 
und der zweiten weiteren Phase beriicksichtigt , ausgeglichen 

jL S ti • 

9. Modellierverfahren nach einem der obigen Anspriiche, 
dadurch gekennzeichnet, dass das Ste- 
fan-Problem eindimensional angesetzt und gelost wird bzw. die 
Stefan-Probleme eindimensional angesetzt und gelost werden 
und dass der Anteil (p3) an Austenit anhand einer nicht line- 
aren Funktion der Lage der Phasengrenze (11, 12) bzw. der 
Phasengrenzen (11, 12) ermittelt wird. 

10. Modellierverfahren nach einem der obigen Anspriiche, 
dadurch gekennzeichnet, dass die Er- 
mittlung, in welchen Konzentrationen (kl bis k4) das mindes- 
tens eine bewegliche Legierungselement beidseits der ersten 
Phasengrenze (11, 12) bzw. beidseits der ersten und beidseits 
der zweiten Phasengrenze (11, 12) vorliegt, anhand der 
Gibbsschen freien Enthalpien (Gl, G2, G3 ) der Phasen erfolgt. 
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11. Modellierverfahren nach einem der obigen Anspriiche, 
dadurch. gekennzeichnet, dass anhand 
der im Anf angszustand (ZA) bereits vorliandenen Phasen und an- 
hand der Gibbsschen freien Enthalpien (Gl, G2, G3 ) der Phasen 
erxnittelt wird, ob sowohl Austenit als auch die erste weitere 
Phase vorliegen bzw. ob zusatzlich zu Austenit und der ersten 
weiteren Phase auch die zweite weitere Phase vorliegt. 

12. Modellierverfahren nach einem der obigen Anspruche, 
dadurch gekennzeichnet, dass das 
Stahlvolximen (1) eine Vielzahl von Volximenelementen (9) um- 
fasst, dass das Stefan-Problem bzw. die Stef an-Probleme nur 
fiir einen Teil der Volumenelemente (9) gelost wird bzw. wer- 
den und dass die lokalen Anteile (pi, p2, p3 ) der Phasen der 
anderen Volumenelemente (9) anhand der lokalen Anteile (pi, 
p2, p3) der Phasen des Tells der Volumenelemente (9) ermit- 
telt werden. 

13 . Modellierverfahren nach Anspruch 12 , 

dadurch gekennzeichnet, dass die 
Warmeleitungsgleichung individuell fiir jedes Volumenelement 
(9) gelost wird. 

14. Modellierverfahren nach einem der Anspriiche 1 bis 13, 
dadurch gekennzeichnet, 

- dass dem Rechner (4) ein Erstzustand (Z) und zumindest eine 
gewunschte EndgroSe (f*) vorgegeben werden, 

- dass das Modellierverfahren nach einem der Anspruche 1 bis 
10 iterativ angewendet wird, 

- dass der Anf angszustand (ZA) der ersten Iteration dem Erst- 
zustand (Z) und der Anf angszustand (ZA) jeder weiteren Ite- 
ration dem unmittelbar zuvor ermittelten Folgezustand (ZF) 
entspricht und 

- dass anhand des nach einer letzten Iteration ermittelten 
Folgezustands (ZF) eine erwartete EndgroSe (f) ermittelt 
und mit der gewunschten EndgroSe (f*) verglichen wird. 
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15. Modellierverfahren nach Anspruch 14, 

dadurch gekennzeichnet, dass es on- 
line und in Echtzeit oder offline ausgefiihrt wird. 

16. Modellierverfahren nach Anspruch 14 oder 15, 
dadurch gekennzeichnet, dass die Ein- 
flussgroBen (W) der Iterationen in ihrer Gesamtheit einem 
Einf lussgrofienverlauf entsprechen, dass der Rechner (4) auf- 
grund des Vergleichs der erwarteten EndgroEe ( f ' ) mit der ge- 
wiinschten EndgroSe (f *) den Einf lussgroiSenverlauf variiert 
und ausgehend vom Erstzustand (Z) das Modellierverfahren nach 
Anspruch 12 erneut ausfuhrt, bis zuiaindest die erwartete End- 
groSe (f) mit der gewiinschten EndgroBe (f*) korrespondiert . 

17. Modellierverfahren nach einem der Anspriiche 1 bis 13, 
dadurch gekennzeichnet, dass es on- 
line \ind in Echtzeit ausgefuhrt wird, dass der Rechner (4) 
anhand einer aus dem Anf angszustand (ZA) ermittelten Anfangs- 
groBe (f) imd einer gewunschten FolgegroBe (f*) die Einfluss- 
grSSe (W) ermittelt und dass der Rechner (4) eine Beeinflus- 
sungseinrichtvmg (2) derart ansteuert, dass das Stahlvolumen 

(1) entsprechend der ermittelten Einf lussgroSe (W) beein- 

flusst wird. 

18 . Datentrager mit einem auf dem Datentrager gespeicherten 
Computerprogramm (6) zur Durchfiihrung eines Modellierverf ah- 
rens nach einem der obigen Anspruche. 

19 . Rechner mit einem Massenspeicher ( 8 ) , in dem ein Compu- 
terprogramm (6) hinterlegt ist, so dass bei Aufruf des Compu- 
terprogramms (6) von dem Rechner ein Modellierverfahren nach 
einem der Anspruche 1 bis 17 ausfiihrbar ist. 

20. Beeinflussiingseinrichtung zum Beeinflussen der Tempera tur 
eines Stahlvolumens (1), insbesondere Kuhlstrecke, 
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dadurch gekennzeichnet, dass sie von 
einem Rechner (4) nach Anspruch 19 gesteuert wird. 

21. Stahl, 

dadurch gekennzeichnet, dass er eine 
Beeinflussungseinrichtung (2) nach Anspruch 20 durchlaufen 
hat, wobei die Einf lussgroSe (W) nach Anspruch 16 oder 17 er- 
mittelt wurde. 
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Rechnergestiitztes Modellierverfahren fur das Verhalten eines 
Stahlvo lumens mit einer Volumenoberf lache 

Ein Rechner (4) ermittelt anhand eines momentanen Anfangszu- 
stands (ZA) eines Stahlvo liamens (1) vLad mindestens einer iiber 
die Volijmenoberf lache auf das Stahlvolumen (1) einwirkenden 
momentanen Einf lussgroSe (W) durch Losen einer Warmeleitungs- 
und einer PhasentJiawandlungsgleichung einen Folgezustand (ZF) 
des Stahlvolumens (1) . Die Zustande (ZA, ZF) umfassen fiir 
mindestens ein Volumenelement (9) des Stahlvolumens (1) eine 
lokale Konzentrationsverteilung (K) eines im Stahl bewegli- 
chen Legierungselements , lokale Anteile (pi, p2 , p3) von mo- 
dellierten Phasen des Stahls vmd eine den lokalen Energiein- 
halt des Stahls beschreibende GroSe (H) . Die Phasen umfassen 
Austenit und eine weitere Phase, in der Kegel Ferrit oder Ze- 
mentit. Im Rahmen der Umwandlixngsgleichung wird ermittelt, 
welche Konzentrationen (kl, k3; k2, k4) des beweglichen Le- 
gierungselements beidseits einer Phasengrenze (11, 12) zwi- . 
schen Austenit und der weiteren Phase vorliegen. Ferner wird 
durch Losen eines Stefan- Problems ermittelt, ob und wie sich 
die Konzentrationsverteilung (K) des beweglichen Legierungs- 
elements im austenitischen Bereich des betrachteten Volumen- 
element s (9) andert und ob und vao. welches AusmaS (5x' , 8x" ) 
sich die Phasengrenze (11, 12) dadurch verschiebt. Anhand der 
Lage der Phasengrenze (11, 12) werden dann die lokalen An- 
teile (pi, p2, p3) der Phasen ermittelt. 
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FIG 4 
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nimm entgegen: Z, f*,- 
Z=Z(p1.p2, p3, H,K) 
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bestimme EinflussgroB enverlauf 



S3 



ZA=Z 

t=to 



S4 



bestimme W(t) 



bestimme anhand von ZA und W durch Losen 
der Warmeleitungsgleictiung und der 
Phasenumwandlungsgieicliung ZF 

ZF=ZF(p1,p2, p3,H,K) 
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t=t+5t 
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ZA=ZF 
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S8 



variiere Einfluss- 
groBenverlauf 



S10 



steuere Kulilstrecl<e entsprecliend dem i 
ermittelten EinflussgroBenverlauf an , 
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nimm entgegen: Z 




Z=Z(p1,p2, p3,H,K) 
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ZA=Z 

t=to 
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nimm entgegen 
Oder ermittle: f* 
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S16 



ermittle f(ZA) 



SI 7 



ermittle anhand von 
f und f*W 



steuere Kuhlstrecke entsprechend an 



bestimme anhand von ZA und W durch Losen 
der Warmeleitungs- und der Phasen- 
umwandlungsgleichung ZF 

ZF=ZF(p1,p2, p3, H,K) 
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FIG 6 
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bilde Gruppen von 
Volumeneiementen 



S40 



FIGS 



S41 



wahle pro Gruppe ein 
Volumenelement aus 



S22 



ermittle fiir ausgewahlte 
Volumenelemente anhand 

H=p1H1(T)+p2H2(T)+p3H3(T) 
lokaleTemperaturT 



ermittle bei dieser Temperatur 

Gibbsche freie Enthalpien 
der Phasen und Konzentrationen 
der Phasen an Phasengrenzen 



S23 
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ermittle Anzahl der zu 
betrachtenden Phasen 
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lose Diffusionsgleichung 




S30 



1 



lose Stefan-Problem 
fiir Austenit und 
weitere Phase 
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lose Stefan-Probleme fiir 
Austen it-Ferrit und fiir 
Austenit-Zementit; beachte 
Kopplung der Stefan-Probleme 



S37 



S29 



bestimme anhand von 6x neue 
Lage der Phasengrenze(n) 
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p3(ZF)=nicht lineare Funktion (Lage) 
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p1 (ZF) = p1 (ZA)-(1 -q). [p3(ZF)-p3(ZA)] 
p2(ZF) = p2(ZA)-q • [p3(ZF)-p3(ZA)] 

i- = 



ubernimm Ergebnis fiir die anderen Volumenelemente 
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lose Warmeleitungsgleichung fiir alle Volumenelemente 
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FIG 9 
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FIG 11 
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FIG 19 
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FIG 18 
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mit L1 = 
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